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INTRODUCTION 

Characterizing  soil  behavior  under  loading  is  a  complex  task.  In 
recent  years  the  finite  element  method  has  been  a  useful  tool  in  analyz¬ 
ing  structures,  structures  on  soil,  and  soil-structures  themselves. 
Emphasis  in  the  past  has  been  placed  on  use  of  nonlinear  plasticity  soil 
models  to  capture  soil  response  more  accurately.  Recent  work  has  fo¬ 
cused  on  effective  stress  analysis  -  the  ability  not  only  to  calculate 
soil  stress  but  also  to  calculate  pore  fluid  pressure. 

An  effective  stress  model  can  be  applied  to  analysis  of  ocean  floor 
soils,  nearshore  and  offshore  structures,  and  seismic  phenomena.  Oscil¬ 
lations  in  loading,  whether  from  wave  action  or  seismic  shaking,  produce 
dynamic  loading  that  can  induce  significant  increases  in  pore  pressure. 
In  turn  this  increase  in  pore  pressure  can  reduce  allowable  soil  capac¬ 
ities  and  at  the  same  time  increase  deformations  from  a  reduction  in 
effective  confining  stress.  Under  extreme  conditions,  flow  slides  and 
liquefaction  occur.  Although  liquefaction  has  been  identified  as  a 
phenomenon  for  20  years,  soil  mechanics  experts  are  just  beginning  to 
understand  the  interaction  of  stress  confinement  and  drainage  path  that 
occurs  in  the  field,  such  as  under  a  foundation  or  around  a  pile.  For 
example,  common  engineering  practice  in  evaluating  seismically  induced 
soil  liquefaction  considers  ground  conditions  as  level  away  from  the 
structure,  and  shear  stresses  from  the  structure  are  not  even  consid¬ 
ered.  Also,  present  design  guidelines  for  pile  foundations  are  based  on 
static  load. 

Recent  earthquakes,  particularly  those  in  Alaska,  Japan,  and  Chile, 
have  emphasized  the  high  damage  threat  the  soil  liquefaction  phenomenon 
poses  to  waterfront  structures.  In  the  1960  Chilean  earthquake 
(magnitude  7)  quaywalls,  sheet  piles,  and  sea  walls  were  damaged  by 


liquefaction  of  loose,  fine,  sandy  soils.  In  the  1964  Alaskan  earth¬ 
quake  (magnitude  8.4)  severe  damage  to  Anchorage,  Cordova,  and  Valdez 
occurred,  including  large-scale  landslides,  as  a  result  of  liquefaction. 
Japanese  earthquakes  (Niigata,  1964,  magnitude  7.3;  Miyagi-Ken-Oki , 
1978,  magnitude  7.4)  caused  severe  waterfront  damage  to  wharfs,  bulk¬ 
heads,  quaywalls,  piers,  and  conventional  structures.  The  majority  of 
the  damage  sustained  in  waterfront  areas  was  from  liquefaction  of  loose, 
cohesionless  sands. 

A  study  conducted  by  the  Office  of  Naval  Research  (Ref  1)  recog¬ 
nized  a  major  liquefaction  hazard  existing  at  West  Coast  Naval  activ¬ 
ities.  A  more  recent  investigation  at  the  Naval  Air  Station  (NAS), 
North  Island,  Calif.  (Ref  2),  concluded  that  liquefaction  under  design 
earthquake  levels  could  result  in  destruction  of  such  critical  struc¬ 
tures  as  aircraft  carrier  berths,  aviation  fuel  tank  farms,  and  under¬ 
ground  utility  service  lines.  Unfortunately,  almost  all  previous 
studies  of  the  liquefaction  problem  have  been  concerned  with  either 
conventional  building  foundations  or  with  analyses  of  dams,  and  proce¬ 
dures  for  analyzing  specialized  Navy  structures  are  not  available. 

The  effective  stress  model  is  of  major  significance  since  the  Navy 
must  locate  in  areas  where  the  water  table  is  high.  Even  if  liquefac¬ 
tion  (a  loss  of  shear  strength  from  a  loss  of  effective  confining 
stress)  does  not  occur,  a  buildup  of  pore  pressure  is  probable  both  in 
sands  and  clays.  This  pore  pressure  buildup  can  be  of  major  signif¬ 
icance  in  structural  behavior. 

The  Navy  has  a  drydock  certification  program  in  progress.  The 
drydocks  are  examples  of  structures  surrounded  by  soil,  often  with  a 
high  water  table.  Soil  loading  causes  wall  deflections  which,  in  turn, 
alter  the  soil  load.  Note  that  even  though  drydocks  are  critical  Navy 
structures,  the  present  analytical  techniques  reflect  the  state-of-the- 
art  of  1950.  Basically,  static  structural  analysis  procedures  are  used 
with  estimates  of  the  soil  pressure.  The  drydock  is  only  one  of  many 
types  of  waterfront  structures  illustrating  the  significance  of  soil- 
structure  interaction. 


The  effective  stress  soil  model  is  a  critical  tool  for  use  on 
waterfront  structures.  Dynamic  analysis  techniques  are  essential  for  a 
realistic  assessment  of  drydock  safety.  Drydocks  are  only  one  applica¬ 
tion;  others  include  quaywalls,  bulkheads,  retaining  walls,  ocean  floor 
structures,  as  well  as  other  structures. 

Effective  stress  techniques  are  also  useful  in  understanding  site 
behavior  because  pore  pressure  builds  up  during  the  cyclic  loading 
process,  increasing  deformations. 

Conventional  finite  element  techniques  presently  utilize  a  total 
stress  approach  characterizing  the  soil  elastically  in  terms  of  shear 
modulus  and  Poisson's  ratio  or  inelastical ly  in  terms  of  a  curved 
stress-strain  fit  typically  from  single  triaxial  compression  and  exten¬ 
sion  tests.  This  approach  does  not  consider  the  pore  fluid  separately 
from  the  soil  skeleton,  but  actual  soil  behavior  depends  heavily  upon 
the  behavior  of  the  fluid.  Initially  upon  loading,  much  of  the  load  is 
carried  by  the  near  incompressible  fluid;  as  drainage  occurs,  load  is 
transferred  to  the  soil  skeleton.  Fluid  flow  controlled  by  soil  perme¬ 
ability  plays  an  important  part  in  computations  of  the  actual  loading 
and  settlements  in  the  soil.  To  neglect  the  pore  fluid,  as  is  presently 
done  in  total  stress  finite  element  programs  used  today,  is  a  gross 
oversimplification  limited  to  dry  soils  and  is  not  well-suited  to  water¬ 
front  construction. 

This  report  summarizes  previous  work,  and  has  been  funded  by  Min¬ 
erals  Management  Service  and  the  Navy.  During  FY81  a  detailed  review  of 
available  material  models  was  made,  and  model  behavior  was  compared  with 
test  data  (Ref  3).  Based  on  the  comparative  study,  it  was  determined 
that  the  Prevost  model  was  applicable  to  both  cohesive  and  cohesionless 
materials  and  that  its  development  was  very  far  along  and  well  into  the 
implementation  stage. 

During  FY82,  modifications  were  accomplished  by  contract  with 
Professor  J.H.  Prevost  to  include  the  effects  of  cyclic  degradation. 
Soil  test  data  were  obtained  by  contract  with  Professor  P.  Lade  of  the 
University  of  California  at  Los  Angeles,  Calif.  (UCLA).  Reference  4 
presents  this  work.  During  FY83,  additional  test  data  were  obtained  by 
contract  with  Dr.  Richard  Ladd,  Woodward- Clyde  Consultants  (Ref  5). 


Test  data  and  model  prediction  were  compared,  and  several  boundary  value 
problems  were  tried  (Ref  6).  During  the  current  fiscal  year,  work 
continued  on  comparison  of  model  prediction  and  test  data  and  also  on 
advanced  boundary  value  problem  formulation  and  solution.  Required 
enhancements  were  made  in  the  program  formulation. 


PREVCST  SOIL  MODEL  FORMULATION 

'In  the  Prevost.soil  model,  soil  is  viewed  as  a  multi-phase  medium 
consisting  of  an  inelastic  porous  skeleton  and  viscous  fluids.  The 
model  (Ref  7  and  8)  is  a  general  analytical  model  that  describes  the 
nonlinear,  anisotropic,  elasto-plastic,  stress  and  strain  dependent,  and 
strength  properties  of  the  skeleton  when  subjected  to  a  three- 
dimensional  loading.  The  following  is  taken  directly  from  Prevost 
(Ref  9). 

The  constitutive  equations  for  the  solid  skeleton  are  written  in 
one  of  the  following  forms: 


a1  small  deformations 

S'  +  o'  div  v  finite  deformations 


(1) 


where:  o'  =  effective  Cauchy  stress  tensor 

v  =  spatial  velocity  of  solid  phase 

£  =  the  rate  of  deformation  tensor  for  the  solid  phase 

(symmetric  part  of  the  spatial  solid  velocity  gradient) 

•  =  material  derivative 

S'  =  Jaumann  derivative,  given  by: 


S' 


=  o'  +  o' 


w  -  w 


(2) 


where  w  is  the  spin  tensor  for  the  solid  phase  (skew- symmetric  part  of 
the  spatial  solid  velocity  gradient).  In  Equation  1,  Cabcd  is  a  tensor 
valued  function  o'  and  the  solid  deformation  gradients. 


For  soil  media,  the  form  of  the  C- tensor  is  given  as  follows, 
(E:P)  <§:JE) 

£  =  s  -  Jr—gTFIp-  <3> 

where:  H'  =  plastic  modulus 

P  and  g  =  symmetric  second-order  tensors,  such  that  P  gives 
the  direction  of  plastic  deformations  and  §  the 
outer  normal  to  the  active  yield  surface 

E  =  fourth-order  tensor  of  elastic  moduli,  assumed 
isotropic,  viz., 

E  .  .  =  A  6  .6.  +  G  (6  6.  .  +  6  A  )  (4) 

abed  ab  be  v  ac  bd  ad  be'  '  J 


where  A  and  G  are  Lame's  constants,  and  6ab  is  the  Kronecker  delta.  The 
yield  function  is  selected  of  the  following  form 


f  =  f  (S  -  a)  :  (S  -  a)  +  C2  (p1  -  p)2  =  k2  (5) 

where:  S  =  a'  -  p*l  *  deviatoric  stress  tensor 

p1  =  (1/3)  tr  o'  =  effective  mean  normal  stress 

a  and  p  =  coordinates  of  the  center  of  the  yield  surface  in  the 
deviatoric  stress  subspace  and  along  the  hydrostatic 
stress  axis,  respectively 

k  =  size  of  the  yield  surface 

C  =  material  parameter  called  the  yield  surface  axis  ratio 


From  Equation  5 


grad  /  =  j£r  =  3  (S  -  a)  ♦  f  -  C2  (p-  -  p)  1 


and 


grad  f 


2  =  6k2  +  6C2  (|  C2  -  1)  (p'  -  p)2 


(6) 


(7) 
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It  is  convenient  to  decompose  P  and  9  into  their  deviatoric  and 
dilational  components,  and  in  the  following 

P  =  P'  +  P"  1  g  =  Q'  +  Q"  1  (8) 

where 

P"  =  |  tr  P  Q"  =  |  tr  9  (9) 

and 

9  =  grad  f/|grad  f  (10) 

The  plastic  potential  is  selected  such  that  the  plastic  rate  of 
deformation  vector  remains  normal  to  the  projection  of  the  yield  surface 
onto  the  deviatoric  stress  subspace 

r  =  s' 

and 

3P"  =  Q"  +  A  tr  (9' )3/tr  (9' )2 

where 

tr  (9'  )2  =  9' : g'  =  g;b  g;b  d2) 

tr(9')3  =  3  det  (91 )  =  lib  Qbc  Qca 

where  A  is  the  material  parameter  which  measures  the  departure  from  an 
associative  plastic  flow  rule.  When  A  =  0,  the  principal  directions  of 
P  and  9  coincide  and  consequently  the  C  tensor  possesses  the  major 
symmetry  and  leads  to  a  symmetric  material  tangent  stiffness.  On  the 
other  hand,  when  A  t  0,  the  principal  directions  of  P  and  9  do 
coincide,  and  C  does  not  possess  the  major  symmetry. 


(11a) 

(lib) 


From  Equations  4,  8,  and  9, 


Q: JE:  P  =  B  tr  P  tr  Q  +  2G  P1  :Q'  (13) 

where  B  =  A  +  2G/3  =  the  elastic  bulk 
G  =  shear  moduli 
For  the  small  deformation  case: 

o'  =  2G  e  +  (B-fS)  lv  1  (14) 


/  2G  Q'  e  +  B  3Q"  e  \ 

(2G  S'  +  8  3P"  i^H'  +  2G  9 '  Q '  +  B  3Q"3P7 

where  =  tr  e;  or,  equivalently,  in  terms  of  deviatoric  and 
dilational  components: 


/  2G  9*  e  +  B  3Q"  \ 

S  -  2G  e  -  2G  9  ^H,  +  2G  g,.g,  +  g  3q„  3pllj  (15a) 

(  2G  9':e  ♦  8  3QH  ev  \ 

p  8  ev  -  B  3P  ^H,  +  2G  g,.g,  +  g  3Qi.3piij  (15b) 


where  e  =  e  -  e  1  =  deviatoric  rate  of  deformation  tensor.  When  O'  =  0, 
2  ~  ~  v~  •* 

3Q"  =  1  and  Equation  15b  simplifies  to 

'  =  <16> 


and  the  plastic  modulus  H'  thus  plays  the  role  of  a  plastic  bulk  mod¬ 
ulus.  Similarly,  when  Q"  =  0,  Q':Q'  =  1  and  Equation  15  yields 


3':!  =  9'  =  S 


(17) 


and  the  plastic  modulus  H'  thus  plays  the  role  of  a  plastic  shear 
modulus. 


When  C  =  0,  the  yield  surface  plots  in  stress  space  as  a  cylinder 
(translated  Von-Mises  yield  surface)  whose  axis  is  parallel  to  the  space 
diagonal.  When  C  t  0,  the  yield  surface  plots  in  stress  space  as  an 
ellipsoidal  surface  of  revolution  whose  axis  is  parallel  to  the  space 
diagonal.  In  order  to  allow  for  the  adjustment  of  the  plastic  hardening 
rule  to  any  kind  of  experimental  data  (for  example,  data  obtained  from 
axial  or  simple  shear  soil  tests),  a  collection  of  nested  yield  surfaces 
is  used.  The  surfaces  are  all  similar;  i.e.,  the  axis  ratio  C  in  Equa¬ 
tion  5  is  the  same  for  all  yield  surfaces.  The  yield  surfaces,  in 
general,  may  translate  and  change  in  size,  but  never  rotate.  The  model, 
therefore,  combines  properties  of  isotropic  and  kinematic  plasticity. 
In  order  to  avoid  overlappings  of  the  surfaces  (which  would  lead  to  a 
non-unique  definition  of  the  constitutive  theory)  the  isotropic/ 
kinematic  hardening  rule  couples  the  simultaneous  deformation/ 
translation  of  all  yield  surfaces.  A  plastic  modulus  H1^  and  a  non- 
associative  parameter  are  associated  with  each  yield  surface,  m. 

In  general,  both  A^  and  H'^  are  allowed  to  take  different  values  at 
different  locations  on  any  given  yield  surface;  i.e.,  both  A^  and 
H1^  are  functions  of  position. 

The  following  rules  have  been  used  for  A^  and 

Cohesive  soils:  A'm^  is  constant  on  each  surface,  and 


H,(m) 


tr 

nTT 


g.  (m) 


(18) 


where  h1^  =  plastic  shear  modulus,  and  [h'^m^  ±  B'^]  =  the  plastic 
bul k  modul i . 

Cohesionless  soils:  Both  A^  and  H'^  are  allowed  to  vary  on 
each  yield  surface; 


A(m)  = 


A^"0 


i=l  if  tr  (g,(m))3  >  0 
i  =2  if  tr  (g,(in))3  <  0 


(19) 
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and 


i  (m)  _  ^  i  (m) 


tr  (S,(m))3 


+  B; 


.  (m) 


1  - 


tr  (Q,(m))3 


(20) 


i=l  if  tr  Q(m)  >  0 
i=2  if  tr  <  0 

where  hi^  and  B^m^  play  the  role  of  plastic  shear  and  bulk  moduli, 
respectively. 

The  yield  surfaces'  initial  positions  and  sizes  reflect  the  past 
stress-strain  history  of  the  soil  skeleton,  and  their  initial  positions, 
in  particular,  are  a  direct  expression  of  the  material's  "memory"  of  its 
past  loading  history.  Because  the  a's  are  not  necessarily  all  equal  to 
zero,  the  yielding  of  the  material  is  anisotropic.  Direction  is  there¬ 
fore  of  importance,  and  the  physical  reference  axes  (x,y,z)  are  fixed 
with  respect  to  the  material  element  and  specified  to  coincide  with  the 
reference  axes  of  consolidation.  For  a  soil  element  whose  anisotropy 
initially  exhibits  rotational  symmetry  about  the  y-axis, 

a  =  a  =  -of  /2,  and  Equation  5  simplifies  to 
x  z  y 

t(°y  -  o;)  -  a]2  +  C2  (p'  -  P)2  -  k2  =  0  (21) 


where  a  is  3a^/2.  The  yield  surfaces  then  plot  as  ellipses  in  the 
axisymmetric  stress  plane  (a^  =  o^)  as  shown  in  Figure  la.  Points  C  and 
E  on  the  outermost  yield  surface  define  the  critical  state  conditions 
(i.e.,  H'  =  0)  for  axial  compression  and  extension  loading  conditions, 
respectively.  It  is  assumed  that  the  slopes  of  the  critical  state  lines 
0C  and  0E  remain  constant  during  yielding. 

The  yielding  surfaces  are  allowed  to  change  in  size  as  well  as  to 
be  translated  by  the  stress  point.  Their  associated  plastic  moduli  are 
also  allowed  to  vary;  in  general,  both  k  and  H'  are  functions  of  the 
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plastic  strain  history.  They  are  conveniently  taken  as  functions  of 
invariant  measures  of  the  amount  of  plastic  volumetric  strains  or  plas¬ 
tic  shear  distortions,  respectively.  Complete  specification  of  the 
model  parameters  requires  the  determination  of: 

(1)  the  initial  positions  and  sizes  of  the  yield  surfaces  together 
with  their  associated  plastic  moduli 

(2)  their  size  or  plastic  modulus  changes  as  loading  proceeds 

(3)  the  elastic  shear  G  and  bulk  B  moduli 


The  yield  surface  is  chosen  as  a  degenerate  yield  surface  of 
size  k^  =  0  which  coincides  with  the  stress  point.  Further,  in  order 
to  get  a  smooth  transition  from  the  elastic  into  the  plastic  regime, 
=  0  and  H'^  =  »,  so  that  the  material  behavior  inside  is 

purely  elastic.  The  dependence  of  the  model  parameters  upon  the  effec¬ 
tive  mean  normal  stress  and  volumetric  strain  are  assumed  of  the  fol¬ 
lowing  form 


respectively,  where 


x  =  B,  G  and  H'(m) 
and 

y  =  «<;>,  p<">  and  k<*> 

where:  n  =  experimental  parameter  (n  =  0.5  for  most  cohesionless 

soils  and  n  =  1  for  most  cohesive  soils) 

pi  =  reference  effective  mean  normal  stress  (i.e. ,  at 
ey  =  0  when  p’  =  p^). 

It  is  assumed  that  when  the  soil  is  in  a  "normally  consolidated"  state, 
the  consolidation  soil  test  results  plot  (1)  as  a  straight  line  parallel 
to  the  projections  of  the  critical  state  lines  in  the  £n(p'/p^)  versus 
ev  diagram  and  (2)  as  a  straight  line  in  the  axial  stress  plane.  The 
parameter  \  is  then  simply  determined  from  the  results  of  Kq  minus 
consolidation  soil  test  results. 


(22b) 


where  the  subscript/superscript  K  refers  to  KQ  minus  loading  conditions. 
The  soil's  anisotropy  originally  develops  during  its  deposition  and 
subsequent  consolidation  which,  in  most  practical  cases,  occurs  under  no 
lateral  deformations.  In  the  following,  the  y-axis  is  vertical  and 
coincides  with  the  direction  of  consolidation,  the  horizontal  xz-plane 
is  thus  a  plane  of  the  material's  isotropy;  and  the  material's  aniso¬ 
tropy  initially  exhibits  rotational  symmetry  about  the  vertical  y-axis 
The  model  parameters  required  to  characterize  the  behavior  of  any  given 
soil  can  then  be  derived  entirely  from  the  results  of  conventional 
monotonic  axial  and  cyclic  strain-controlled  simple  shear  soil  texts. 
The  above  material  in  this  section  is  based  on  Reference  9. 


DETERMINATION  OF  MODEL  PARAMETERS 

Compressive  stresses  and  strains  are  considered  positive  and  all 
stresses  are  effective  stresses  unless  otherwise  specified.  For  a 
material  which  initially  exhibits  cross-anisotropy  about  the  vertical 
y-axis,  the  initial  position  in  stress  space  of  the  yield  surfaces  are 
defined  by  the  sole  determination  of  the  two  parameters  and  0^ 

and  Equation  5  simplified  to 

[q  -  a(m)]2  ♦  C2  [p'  -  p(m)J2  -  [k(lB)]2  =  0  (23) 


for  axial  loading  conditions  (i.e.  ,  o'  =  o'  and  x  =  x  =  x  =  0), 

A  L  AV  V  L  iA 

where  q  =  (o'  -  o').  The  yield  surfaces  then  plot  as  circles  in  the  q 

y  * 

versus  Cp '  plane  (referred  to  as  the  axial  stress  plane  hereafter)  as 
shown  in  Figure  lb.  When  the  stress  point  reaches  the  yield  sur¬ 


face  /w, 

q  =  a(m)  ♦  k{m)  sin  0 


P*  =  pO  +  Jocose 


(24a) 


(24b) 


where  0  is  defined  in  Figure  lb  and  Equation  15  simplifies  to: 


e  _  1  1  sin  0  (sin  0  *  Cy  cos  9) 

q  2G  H'^  sin2  0  +  §  C2  cos2  0 


(25a) 


.  I£ 


cos  0  +  sin  0  4—  (sin  0  +  Cy  cos  0) 


p* 

H1  v  y  L 

sin  0  +  |  C  cos  0 

in  which 

MU 


Figure  1(b).  Prevost  soil  model  -  Yield 
surfaces  in  q-cp  plane. 

Monotonic  Drained  Axial  Compression  and  Extension  Soil  Test  Results 


Let  0r  and  0C  denote  the  values  of  0  when  the  stress  point  reaches 
the  yield  surface  fK  1  in  axial  compression  and  extension  loading  condi¬ 
tions,  respectively.  Combining  Equations  24  through  28,  one  finds  that: 


tan  (0r  +  0p)  =  -~£K  - 

L  b  1  -  R^ 


tan  0, 


3y  fJXV  A  (m) 

JYcWJ  Ac 


(30a) 


j_  =  3_  l/vl  A  <•> 

in  8£  2C  [  ¥e\ye/  E 


(30b) 


in  which 


p£  -  p£  exp  A 

[(4  •  4)J 

qc  -  qE  exp  A 

[(4-4)] 

/„•  \n  : 

l_  _  /fc_\  e  _  JL_ 
XC  \  P1  /  q  2G  1 

1  _  /Pc\nIi  1 

*C  \  »i  /  »'  B1 


and  similarly  for  and  y^,  where  the  subscript  C  and  E  refer  to  axial 
compression  and  extension  loading  conditions,  respectively.  Further, 


,.<•)  -  x  i|11/ll,Vc,i cos  9C  \ 

c  c  <\,in2  ec  -  |  c2  COS2  ej 


and,  similarly,  for  \ 

The  smooth  experimental  stress-strain  curves  obtained  in  axial 
tests  are  approximated  by  linear  segments  along  which  the  tangent  (or 
secant)  modulus  is  constant.  The  degree  of  accuracy  achieved  by  such  as 
representative  of  the  experimental  curves  is  directly  dependent  upon  the 
number  of  linear  segments  used.  The  model  parameters  associated  with 
the  yield  surface  are  determined  by  the  condition  that  the  slopes 
q/e  are  to  be  the  same  in  axial  compression  and  extension  tests  when  the 
stress  point  has  reached  the  yield  surface  f^.  The  corresponding 
values  of  0r  and  0C  are  determined  by  combining  Equations  29  and  30  once 
a  rule  has  been  adopted  for  A'  .  Once  6~  and  0^  have  been  determined, 
the  model  parameters  associated  with  f'ln'  are  simply  obtained  from 
Equations  24  and  25. 


*  w‘-‘ 
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Monotonic  Undrained  Axial  Compression  and  Extension  Soil  Test  Results 


In  undrained  tests,  ew  =  0,  and  (from  Equation  33)  yr  =  yc  =  -B,  in 
that  case.  The  model  parameters  associated  with  the  yield  surface  1 
are  again  determined  by  the  condition  that  the  slopes  q/e  are  to  be  the 
same  in  axial  compression  and  extension  tests  when  the  stress  point  has 
reached  the  yield  surface  f^.  As  previously,  the  corresponding  values 
and  0£  are  determined  from  Equations  29  and  30,  in  which 


R  =  C 


P  r 


“  ^E 


(35) 


Knowing  6r  and  0C,  the  model  parameters  associated  with  are 

L  t  C  E 

computed  from  Equations  24  and  25,  in  which  e— =  e— =  0. 

Simple  Shear  Soil  Test  Results 

In  simple  shear  soil  tests,  e  =  e  =  e  =  0.  The  necessary  alge- 

x  y  z 

bra  for  the  determination  of  the  model  parameters  is  considerably  sim¬ 
plified  in  that  case  if  the  elastic  contributions  to  the  normal  strains 
are  neglected.  Equation  15  then  yields  a'x  =  o^,  and  it  follows  that: 


(m) 

or 

p(m) 

k(m) 


= 

=  -  f  (a;  -  a;) 

=  ^l^y 


(36) 

(37) 

(38) 

(39) 


The  model  parameters  associated  with  the  yield  surface  '  are 
then  simply  determined  from  the  above  equations  and  an  incremented 
linear  representation  of  the  shear  stress-strain  curves  obtained  from  a 
simple  shear  test.  Note  that  the  sole  use  of  simple  shear  test  results 
does  not  allow  the  determination  of  the  parameters  B'^  and  On 
the  other  hand,  it  is  apparent  from  Equations  36  and  39  that  the  degra¬ 
dation  of  the  mechanical  properties  of  the  material  under  cyclic  shear 
loading  conditions;  i.e.,  k^  (e)  and  h1^  (e)  with 


*/(iM 


k  =  £  -  3  (tr  £)  l 


where  the  integration  is  carried  along  the  strain  path,  are  most  con¬ 
veniently  determined  from  the  results  of  cyclic  strain-controlled  simple 
shear  tests 


(e  = 


fh 


in  that  case). 


This  is  explained  and  further  discussed  in  Reference  8.  The  above 
material  is  based  directly  on  Prevost  (Ref  9). 


SCALING  DATA  FOR  DEPTH 


Based  on  limited  data,  it  has  been  found  reasonable  to  scale  by  the 
ratio  of  average  confining  stress  levels.  The  elastic  shear  modulus, 
the  initial  elasto-plastic  shear  modulus,  and  the  elastic  bulk  modulus 
are  scaled  by  the  square  root  of  the  stress  ratio.  The  bulk  exponent  is 
not  scaled.  Initial  stress  components  are  scaled  by  the  ratio.  The 
softening  parameters  --  delta,  delta  ultimate,  and  the  yield  surface 
axis  ratio  --  are  scaled  by  the  stress  ratio.  The  size  of  the  yield 
surface  is  scaled  by  the  stress  ratio.  The  elasto-plastic  shear  and 
plastic  bulk  moduli  are  scaled  by  the  square  root  of  the  stress  ratio; 
the  degree  of  nonassociativity  is  not  scaled. 
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Anisotropic  consolidation  can  also  be  treated  by  shifting  the 
ellipses.  It  is  assumed  that  isotropic  data  are  used  and  the  value 
is  maintained  constant,  as  follows: 


=  (1  +  2  ko)  av 


F  (o' ,  n)  =  o  (49) 

where  is  from  the  isotropic  test  and  K  is  the  coefficient  of  lateral 
earth  pressure  at  rest. 

It  is  suggested  that  data  be  first  scaled  to  the  correct  and 
then  shifted  to  anisotropic  conditions. 
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DYNAFLOW  FINITE  ELEMENT  CODE 


The  DYNAFLOW  Program  (Ref  10)  is  a  finite  element  program  intended 
for  the  static  and  transient  analysis  of  linear  and  nonlinear  two-  and 
three-dimensional  problems.  The  analysis  capabilities  include  the 
following: 

Static  -  Nonlinear  elliptic  boundary  value  problems  with  two 
degrees-of- freedom  per  two-dimensional  node  and  three  degrees-of- 
freedom  per  three-dimensional  node. 

Diffusion  -  Nonlinear  parabolic  boundary  value  problems  with  one 
additional  degree-of- freedom  for  the  phase  porous  fluid  pore 
pressure. 

Dynamic  -  Nonlinear  hyperbolic  boundary  value  problems  with  two 
additional  degrees-of- freedom  for  the  phase  porous  fluid. 

The  program  incorporates  the  full  Prevost  model  for  nonlinear 
effective  stress  analysis.  The  material  model  produces  a  nonsymmetric 
stiffness  matrix  Solution  techniques  employed  allow  for  implicit, 
explicit,  and  implicit/explicit  matrix  solution. 

For  static  analysis,  an  incremental  predictor-corrector  load-step 
procedure  is  utilized.  The  term  a  =  1  is  utilized  and  associates  the 
midpoint  of  the  load  (time)  interval  to  the  predictor  phase  and  achieves 
second-order  accuracy. 

For  parabolic  analysis  diffusion  problems,  a  Newton-Raphson  itera¬ 
tive  procedure  is  utilized.  Unconditional  stability  is  achieved  for 
or  >  1/2,  and  a  value  of  a  =  1  is  recommended  to  maximize  high  frequency 
numerical  dissipation.  For  explicit  elements,  a  time-step  restriction 
occurs: 


where  A  is  equal  to  the  largest  eigenvalue  associated  with  the  problem. 


For  hyperbolic  analysis,  an  implicit-explicit  predictor- 
multicorrector  algorithm  is  used.  Stability  is  achieved  when  a  >  1/2. 
Implicit  elements  are  unconditionally  stable  if  0  >  ot/2.  Explicit 
elements  have  a  time-step  restriction: 

Q  <  2  (1  -  £)/(a  +  1/2) 

where  Q  is  uuAt,  and  £  is  the  modal  viscous  damping.  When  a  =  3/2  and 
0=1,  the  dissipative  properties  of  the  algorithm  damp  all  dynamic 
transient  phenomena  and  permit  static  solution. 

Solution  of  the  porous  media  problem  requires  defining  boundary 
conditions  for  the  fluid  as  would  occur  in  a  conventional  seepage  anal¬ 
ysis.  This  is  accomplished  by  defining  the  third  and  fourth  degrees  of 
freedom,  the  x-  and  y-fluid  velocities. 

A  standard  4-node  element  is  used.  To  treat  the  initial  pore  fluid 
incompressibility,  it  is  suggested  that  one-point  Gaussian  quadrature  be 
used  for  the  volumetric  stiffness  numerical  integration  and  that  two-by- 
two  Gaussian  quadrature  be  used  for  the  deviatoric  stiffness  numerical 
integration.  Use  of  elastic-plastic  equations  of  the  above  type  in 
analysis  of  boundary  value  problems  requires  that  an  efficient, 
"sturdy,"  and  accurate  numerical  integration  procedure  for  the  plas¬ 
ticity  equations  at  the  stress  point  level  be  available.  Substantial 
efforts  have  thus  been  devoted  to  designing  a  computational  procedure 
with  the  best  balance  of  accuracy  and  computation  speed.  The  integra¬ 
tion  algorithm  presently  used  in  the  stress  routine  is  a  generalization 
of  the  conventional  radial  return  technique. 

MODEL  EVALUATION 

Cohesive  Soils 

The  soil  data  used  in  this  section  is  part  of  the  data  which  were 
collected  by  the  organizing  committee  of  the  NSF/NSERC  North  American 
Workshop  on  plasticity  and  generalized  stress-strain  applications  in 
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soil  engineering  held  May  28-30,  1980  at  McGill  University,  Montreal, 
Canada  and  reported  by  Prevost  (Ref  9).  Laboratory  axial  test  data  on  a 
laboratory-prepared  Kaolinite  clay  had  then  been  provided  to  Prevost. 
Predictions  about  the  constitutive  behavior  of  the  soil  subjected  to 
loading  stress  paths  not  identified  in  the  data  had  been  requested  by 
the  organizing  committee.  This  section  describes  the  test  results  by 
Prevost  and  their  analysis  and  compares  the  model  predictions  with 
observed  behavior  in  the  tests. 

The  experimental  tests  had  been  conducted  on  cylindrical  samples  in 
a  torsional  shear  testing  device.  All  samples  had  first  been  Kq- 
consolidated  with  a  cell  pressure  of  58  psi  and  a  backpressure  of 
18  psi;  in  other  words,  the  excess  axial  load  necessary  for  Kq- 
consol idation  was  then  released.  All  the  tests  were  stress-controlled 
and  performed  under  constant  volume  conditions  (i.e.,  undrained). 

Figure  2  shows  in  dashed  lines  the  experimental  results  obtained  in 
conventional  undrained  monotonic  axial  compression/extension  soil  tests 
and  in  solid  lines  the  design  curves  used  to  determine  the  model  param¬ 
eters  for  that  clay.  Note  that  some  data  close  to  failure  have  been 
ignored  in  selecting  the  design  curves  because  they  are  not  consistent 
with  the  rest  of  the  data.  This  inconsistency  may  be  due  to  experi¬ 
mental  difficulties  in  capturing  failure  states  in  stress-controlled 
testing  devices. 

Figure  3  shows  model  predictions  for  a  shear  test  in  which  the 
major  principal  stress  is  inclined  at  6  =  15  degrees,  relative  to  the 
vertical  axis  of  the  soil  specimen  (Ref  9).  Figure  3  also  shows  a 
comparison  between  predicted  and  observed  behavior  of  the  soil  in  these 
tests.  Note  that  all  the  model  predictions  agree  well  with  the  experi¬ 
mental  test  results  (Ref  9). 

Cohesionless  Soils 


The  model  developed  by  Prevost  is  able  to  represent  drained  or 
undrained  test  data.  Figures  4  and  5  show  comparison  of  drained  test 
data  for  the  Silica  sand  and  Banding  sand.  Figure  6  shows  a  comparison 
of  undrained  test  data  based  on  the  undrained  parameters  of  the  Banding 
sand. 
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SHEAR  STRESS  q  =  (o'  -  o'  )  (ka/cm2)  SHEAR  STRESS 


Figure  5.  Banding  sand,  drained  compression. 


Stress  difference  -  mean  stress. 

Figure  6.  Banding  sand,  undrained  compression. 


Typical  liquefaction  behavior  is  shown  in  Figure  7  where  cyclic 
degradation  and  pore  pressure  buildup  can  be  seen.  Typical  test  results 
and  model  comparison  are  shown  in  Table  1  for  the  Banding  sand. 


Table  1.  Results  and  Comparisons  of  Cyclic  Tests 
With  Undrained  Banding  Sand 


Tests 

Amp  1 i tude 

No.  of  Cycles  to 
Liquefaction 

Measured 

Predicted 

Stress-control  led 

Cyclic  compression 

-1.20  kg/cm2 

4 

2.5 

Cyclic  extension 

+1.20  kg/cm2 

2 

2 

Strain-controlled 

Cyclic  compression 

-0.0013 

=20 

20 

Cyclic  extension 

+0.0013 

=24 
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SOU  COLUMN  STUDY 

A  soil  column  approximately  12  meters  in  depth  was  studied; 
Figure  8  shows  the  mesh  simulating  infinite  layers.  Since  the  mesh  is 
only  one  element  wide,  wave  reflections  are  not  possible.  Vertical 
displacements  are  restrained;  horizontal  forces  are  applied  to  equili¬ 
brate  gravity  effects;  and  horizontal  fluid  flow  is  restrained.  A 
sinusoidal  excitation  of  0.102g  was  applied  to  the  base.  Figure  9  shows 
the  acceleration  response,  illustrating  soil  attenuation  of  the  induced 
base  motion  which  would  be  expected  for  loose  soil.  Figure  10  shows  the 
variation  of  vertical  effective  stress  with  time  at  four  depths.  The 
stress  decreases  with  increasing  pore  pressure.  This  shows  an  increase 
in  vertical  effective  stress  at  the  point  in  time  that  drainage  begins. 
Figure  11  shows  the  pore  pressure  and  Figure  12  the  shear  stress  at  the 
four  depths.  The  expected  buildup  of  pore  pressure  and  drainage  of 
near-surface  layers  are  illustrated  in  the  results  of  this  case  study. 


Figure  9.  Acceleration  response. 
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(c)  Element  3.  (d)  Base,  Element  1 


Figure  12.  Shear  stress. 


SOIL- STRUCTURE  PROBLEM 


Static  Analysis 

To  illustrate  a  typical  soil-structure  problem,  beam  elements  were 
used  to  represent  a  structural  frame  on  a  soil  foundation.  Figure  13 
shows  the  mesh.  The  structure  was  embedded  in  the  soil  resting  on  the 
top  of  the  fourth  soil  layer.  Figure  14  illustrates  the  coupling  of  the 
beam  elements  to  the  soil  field  using  contact  elements.  Contact  ele¬ 
ments  were  used  at  the  base  to  permit  change  from  the  third  degree-of- 
freedom  in  the  quadrilateral  elements  (pore  pressure)  to  the  third 
degree-of- freedom  in  the  beam  (rotation).  Vertical  and  horizontal 
contact  elements  were  used  along  the  side  of  structure  --  vertical  to 
allow  for  frictional  forces  and  horizontal  to  permit  contact/separation 
with  deformation.  Figure  15  shows  the  deformed  mesh  at  steps  5  and  15 
of  loading.  Figure  16  shows  the  vector  velocity  at  steps  5  and  15. 
Figure  17  shows  the  stress  path  in  the  element  beneath  the  corner  of  the 
structure.  Note  that  as  the  load  is  applied,  the  stress  builds  up 
gradually;  and,  as  shown  in  Figure  18,  most  of  the  load  is  carried  by 
pore  pressure.  As  the  pore  pressure  dissipates  (Figure  18a),  the  effec¬ 
tive  stress  in  the  soil  increases,  and  deformation  also  increases 
(Figure  18b).  The  high  loading  produces  localized  yielding.  This  is 
evident  in  the  deformed  mesh  plots  and  vector  velocity  plots  (Figures  15 
and  16),  which  show  (to  an  exaggerated  scale)  the  failure  of  the  soil 
near  the  structure.  Figure  19  shows  the  pore  pressure  contours.  Pore 
pressure  is  at  maximum  near  the  lower  left  corner  of  the  structure 
because  drainage  is  permitted  only  at  the  surface. 


Displacement  (cm)  Pore  Pressure  (kg/cm 

-3.0  -2.0  -1.0  0.  1.0  0.  0.5 


(a)  Pore  pressure  history  (Element  68). 


(b)  Displacement  history  (Node  85). 

Figure  18.  Two-dimensional  soil  field  with  building. 
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Figure  19.  Pore  pressure  contours  -  Prevost  soil  model. 


Dynamic  Analysis 


It  was  of  interest  to  analyze  a  two-dimensional  model  of  a  struc¬ 
ture  on  soil  subject  to  horizontal  shaking  motion  (see  Figure  20).  A 
full  mesh  is  utilized  to  capture  asymmetric  effects,  such  as  building 
rocking.  The  model  is  composed  of  120  quadrilateral  elements,  43  beam 
elements,  and  33  nodal  contact  elements.  Five  soil  layers  were  used, 
each  composed  of  a  14-surface  model  scaled  to  each  layer's  center.  The 
Prevost  model  incorporates  an  initial  stress  as  a  model  parameter; 
gra  ity  was  applied  from  time  equal  to  0  to  offset  the  initial  material 
stress.  The  building  weight  was  applied  slowly  in  five  steps,  and  pore 
pressures  from  this  load  were  allowed  to  dissipate  to  produce  an  equi¬ 
librium  static  condition  as  the  starting  point  for  dynamic  excitation. 
This  solution  was  performed  using  the  hyperbolic  option  with  integration 
parameters  a  =  1.5  and  0  =  1.0  to  produce  a  highly  damped  "static" 
condition.  Thus,  the  first  five  steps  were  a  static  analysis.  The 
solution  utilized  four  degrees- of- freedom  per  node.  The  vertical  side 
boundaries  were  free  boundaries  for  translation,  having  static  gravity 
nodal  forces  loads  and  fluid  nodal  forces  applied  horizontally  to  main¬ 
tain  equilibrium  and  having  nodal  dampers  to  simulate  nonreflecting 
boundaries.  The  water  table  was  at  the  ground  surface,  where  drainage 
was  permitted.  Figure  21  shows  the  deformed  mesh  at  the  end  of  the 
static  solution  (fifth  step). 


Figure  20.  Undeformed  mesh. 
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(a)  End  of  static  solution,  step  5. 


(c)  Step  110. 


Figure  21.  Deformed  mesh  plots. 
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At  the  sixth  step,  the  time  increment  was  changed  to  a  short  dy¬ 
namic  increment  (0.025  seconds),  and  the  integration  parameters  a  and  p 
were  changed  to  0.65  and  0.33,  respectively,  for  dynamic  damped  re¬ 
sponse.  The  loading  chosen  for  this  case  was  a  sine  wave  horizontally 
applied  at  the  base  as  a  displacement  function.  Since  the  solution  was 
a  continuing  solution  at  the  sixth  step,  initial  conditions  corre¬ 
sponding  to  an  initial  base  velocity  could  not  be  input  during  the 
solution.  One  approximate  technique  to  eliminate  the  effect  of  the 
initial  condition  is  to  merge  the  desired  sinusoidal  excitation  loading 
with  a  cubic  function  such  that  derivatives  of  the  function  are  ini¬ 
tially  0.0.  To  eliminate  the  fluid  pressure  buildup  at  the  boundaries, 
the  boundaries  were  defined  to  permit  flow  across  them.  Fluid  forces 
were  computed  based  on  the  depth  to  the  center  of  each  element  and 
applied  to  the  nodal  points  acting  on  the  third  degree-of- freedom,  as 
follows: 

Fw  =  <y  z  n)  x  element  height 

where:  Fw  =  total  force  applied  (one-half  to  each  of  two  boundary  nodes' 
third  degree- of- freedom  with  sign  corresponding  to 
direction) 

Yw  =  unit  weight  of  fluid 
z  =  depth  to  center  of  element 
n  =  porosity 

The  forces  active  to  restrain  the  vertical  boundary  from  the  soil  weight 
were  computed  and  applied  to  the  first  degree-of- freedom  as  follows: 

Fg  =  [ys  z  (1  -  n))  x  element  height 

where:  F  =  total  force  applied  one-half  to  each  of  the  two  boundary 
nodes'  first  degree-of-freedom  with  sign  corresponding 
to  direction 


Y,  =  unit  weight  of  solids 


Figures  21b  and  21c  presents  the  deformed  mesh,  showing  the  soil  field 
and  structure  displacement  to  an  exaggerated  scale.  Figure  22  shows 
velocity  vector  plots.  Figure  23  shows  the  out-of-balance  forces,  which 
indicate  the  stability  of  the  solution.  The  solution  reaches  a  stable 
level  of  force  imbalance  oscillating  around  400  kg.  Computational 
stability  was  not  a  problem  for  the  time  step  selected  and  probably 
could  be  continued  several  fold  without  increased  error.  Figures  24 
and  25  show  acceleration  and  velocity  in  the  free  field,  near  the  struc¬ 
ture  and  at  the  base  (input  motion)  (Figure  20  shows  nodal  number  loca¬ 
tions).  The  results  show  minor  amplification  of  motion.  Figure  26  shows 
the  mean  effective  stress  after  the  static  solution,  step  5.  Figure  27 
shows  the  mean  effective  stress  at  points  in  the  solution;  Figure  28 
shows  the  pore  pressure  at  the  same  points  in  time.  Figure  29  gives 
pore  pressure  and  mean  stress  histories  at  two  locations  (see 
Figure  20). 

Results  of  the  analysis  show  increases  in  pore  pressure  around  the 
sides  of  the  structure  and  beneath  the  structure  with  horizontal 
shaking.  The  level  of  increase  is  consistent  with  the  O.lg  excitation. 

The  problem  was  repeated  with  a  base  acceleration  of  0.25g.  In 
this  case,  attentuation  of  ground  motion  resulted.  Pore  pressures 
increased  over  those  of  the  previous  run. 


PILE  AND  PIER  FOUNDATIONS 
Load  Capacity 

Piles  and  piers  are  used  to  transmit  foundation  loads  to  strata  of 
adequate  bearing  capacity  and  to  eliminate  settlements.  Piles  also  are 
often  used  to  resist  lateral  loads  or  uplift  forces.  Pile  use  often 
determines  the  type  of  pile  chosen:  timber  piles  are  well-suited  for 
use  as  friction  piles  in  granular  materials;  steel  piles  are  best  for 
endbearing  on  rock;  concrete  cast-in-place  piles  are  economical  to  use; 
and  concrete-filled  steel  piles  offer  high  bending  resistance.  Depend¬ 
ing  on  the  type  of  pile  used,  piles  are  either  driven  or  cast-in-place. 
Deep  foundations  utilize  piers  that  resemble  large  cast-in-place  piles 
to  transmit  the  loading. 


400  - 


Building 


(d)  Step  110. 

Figure  28.  Pore  pressure  beneath  bulding,  kg/cm2. 


Note:  Gravity  effects  removed  (minus  indicates 
below  original  condition). 
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(a)  Element  3. 


(b)  Element  12. 


Figure  29.  Pore  pressure  time  histories 


This  discussion  will  consider  only  straight  piles  driven  into 
homogeneous  deposits  of  cohesionless  materials.  Experience  has  shown 
that  when  piles  are  driven  into  sand,  the  soil  near  the  pile  is  com¬ 
pacted  to  a  distance  of  a  few  pile  diameters.  In  a  homogeneous  sand  the 
point  resistance  and  average  skin  friction  increase  with  depth  of  pene¬ 
tration  up  to  a  critical  depth.  Beyond  the  critical  depth,  the  point 
resistance  and  skin  friction  remain  almost  constant;  generally,  this  is 
caused  by  soil  compressibility,  crushing,  and  arching.  The  empirical 
approach  to  prediction  of  pile  behavior  has  proven  more  satisfactory 
than  an  analytical  bearing  capacity  approach. 

The  load  transfer  mechanism  between  the  pile  and  the  surrounding 
soil  governs  the  behavior  of  the  pile.  The  design  of  a  pile  requires 
determination  of  the  proportion  of  load  transferred  to  the  soil  by 
adhesion  and  friction  between  the  pile  and  the  soil  and  that  transferred 
by  end  bearing.  This  is  influenced  by  the  flexibility  of  the  pile,  the 
stiffness  of  the  soil,  and  the  nature  of  the  transfer  mechanisms  between 
pile  surface  and  soil. 

Through  use  of  instrumented  field  and  model  studies,  the  design  of 
piles  and  piers  has  been  better  understood.  Figure  30  illustrates  use 
of  pile-driving  resistance  formulas  to  estimate  individual  allowable 
pile  loads.  Figure  31  illustrates  the  calculation  of  ultimate  load 
capacity  of  piles  in  cohesive  soils,  and  Figures  32  and  33  illustrate 
the  calculation  of  load  capacity  of  piles  for  cohesionless  soils. 
However,  these  conventional  ultimate  design  approaches  assume  the  simul¬ 
taneous  and  full  mobilization  of  pile  shear  resistance  and  base  bearing, 
which  is  not  well-founded.  Studies  have  shown  that  movement  of  a  pile 
must  be  present  to  mobilize  its  load-carrying  capacity. 

Reese  and  O'Neill  (Ref  11)  show  the  division  of  load  between  pier 
and  base  for  a  pier  on  stiff  clay  (Figure  34).  The  first  two  load 
increments  show  that  almost  all  the  load  is  carried  by  friction;  as  the 
load  increases,  more  is  carried  by  end  bearing.  Above  80  tons,  any 
additional  load  is  taken  by  end  bearing.  Note  how  curve  C  simply  trans¬ 
lates  to  form  curve  D,  indicating  additional  load  is  taken  by  end 
bearing.  With  the  information  in  Figure  34,  a  typical  load  transfer 
relationship  can  be  obtained  showing  side  friction  (Figure  35).  From 
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the  slope  in  Figure  34,  it  is  evident  that  frictional  effects  are  great¬ 
est  in  the  middle  (depth)  of  the  pile;  and  a  reduction  in  the  rate  of 
load  transfer  occurs  in  the  lower  part,  particularly  just  above  the 
base.  The  distribution  of  frictional  forces  depends  upon  the  soil  type. 
Modifications  to  ultimate  strength  formulations  have  been  suggested  by 
Reese  and  O'Neill  (Ref  11).  Factors  include  an  effective  depth  concept 
in  lateral  earth  pressure  calculation  for  granular  soils  and  a  friction 
reduction  factor  for  ultimate  resistance  in  cohesive  soils. 

A  theoretical  solution  of  the  load  transfer  phenomenon  can  provide 
considerable  insight.  For  instance,  a  soil-structure  problem  would  be 
one  in  which  the  nature  and  mechanism  of  load  transfer  is  of  consid¬ 
erable  importance.  Solution  of  this  problem  must  focus  on  the  stress- 
strain  characteristics  of  the  soil  and  the  behavior  of  the  interface 
between  pile  and  soil.  Use  of  a  linear-elastic  constitutive  model  would 
be  a  gross  simplification  of  the  real  material  properties.  The  Prevost 
soil  model  offers  an  opportunity  to  explore  nonlinear  relationships. 

Several  approaches  have  been  used  to  model  the  contact  problem. 
Peterson  (Ref  13)  treats  the  two  contacting  surfaces  as  separate  and 
distinct  and  joins  them  mathematically  by  use  of  Lagrange  multipliers. 
However,  this  does  not  allow  pre-slip  deformation.  Herrmann  (Ref  14, 
15,  and  16)  defines  three  behavior  modes:  nonslip,  slip,  and  separa¬ 
tion.  A  compatibility  model  combining  compatibility  and  equilibrium  is 
used.  The  compatibility  model  involves  linking  the  two  surfaces  with 
fictitious  bond  springs,  and  the  frictional  forces  are  applied  as  sur¬ 
face  tractions.  This  allows  a  pre-slip  deformation  that  would  not 
otherwise  be  computed. 

Another  approach  simply  links  two  nodes,  initially  overlapping  with 
nonlinear  springs.  The  procedure  has  the  advantage  of  simplicity  of 
operation  but,  unfortunately,  may  produce  undesirable  numerical  charac¬ 
teristics  when  distinct  rapid  changes  occur  in  stiffness. 

Prevost  has  formulated  a  contact  element  in  the  DYNAFLOW  code.  The 
contact  element  may  be  used  to  impose  inequality  constraints  between 
nodes.  Either  perfect  friction  (i.e.,  "stick")  or  frictionless  ( i . e .  , 
"slip")  conditions  may  be  achieved. 
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For  single- at  tine  hamm-  r 


For  double- ac  tine; 
oitferenti  i  hammer 


juse  when  driven  QalF  -  .  5 e  when  driven 

weights  are  smaller  5  OJ  weights  are  smaller 

than  striking  than  striking 

weights.  weights. 


2  W  H  [  2  E:  i 

IrfJ  =  - — —  {Use  when  driven  QQf - — '  (Use  when  driven 

S  O  /  weights  are  larger  $  +  Q  / £-  weights  are  larger 

than  striking  '  /Ve  than  striking 


allowable  pile  load  in  pounds 

weight  of  striking  parts  of  hammer  in  pounds 

the  effective  height  ot  fall  in  feet 

the  actual  energy  delivered  by  hammer  per  blow  in  foot-pounds, 
average  net  penetration  in  inches  per  blow  for  the  last  6  in.  of  driving 

driven  weights  7  Note  Ratio  of  driven  weights  to  striking  weights  should  not 

weights  of  striking  parts J  exceed  3. 


Modi :u  ations  ot  basic  pile  driving  tormuius 


A.  For  piles  driven  to  and  seated  in  rock  as  high  capacity  end-bearing  piles: 

Drive  to  refusal  (approximately  4  to  5  blows  for  the  last  quarter  inch  of  driving). 
Redrive  open  end  pipe  piles  repeatedly  until  resistance  for  refusal  is  reached 
within  1  in.  of  additional  penetration. 


B  Piles  driven  through  stiff  compressible  materials  unsuitable  for  pile  bearing  to  an  under¬ 
lying  bearing  stratum: 

Add  blows  attained  before  reaching  bearing  stratum  to  required  blows  attained  in 
bearing  stratum  (see  example). 

grade — -"i  „ 

- - ||  "  Example  :  required  load  capacity  of  pile  «2 5  tons 

F> LL  _ II _ __  hammer  energy  E  =  15,000  ft-lb 


COMPRESSIBLE 

„  STRATUM/ 

W///7/ 


SEARING 
STRA  TUM 


-18  blow*  /  ft 


4?  fclovt  /If 


Penetration(s)  as  per  basic  formuia  =  or  2  blows  per 
inch  (24  blows/ft). 

Required  blows  for  pile  24  -f  18  -  42  blows/ft. 


C  Piles  driven  into  limited  thin  bearing  stratum,  drive  to  predetermined  tip  elevation. 

Determine  allowable  load  by  load  test. 


STRATUM 
UNSUI TABLE 
FOR  BEARING 


5  -  STRATUM 


/STIFF  ClAV  STRATUM  iMCOMPttt  $Sl«LC 

^BUT  UNSUI  TABLE  FOR  POINT  BEAR  INC 

'  /  /  /  /  / 


Figure  30.  Application  of  pile-driving  resistance  formulas  (Ref  11). 
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RECOMMENDED  VALUES  Of  ADHESION 


POSIT IOP  OF  GROUND 
9 A  TER  HAS  90  EFFECT 
09  ULTIMATE  LOAD 
CAPACITY  U9LESS 
C0HESI09  IS  CHA9&E0 

<?„//  IS  APPLIED  LOAD 
09L T . PILE  9EI6HT  IS 
BALA9CED  BT  WEIGHT 
OF  OVERBURDE 9  A90 
IS  90T  C09S I OERED ■ 


ULTIMATE  LOAD  CAPACITY 


PILE 

TYPE 

CONSISTENCY 
OF  SOIL 

COHESION.  C 
PSF 

ADHESION.  CA 
p5F 

VERY  SOFT 

0  -  is 0 

0  -  250 

TIMBER 

SOFT 

ISO  ■  BOO 

250  -  480 

AND 

CONCRETE 

NED.  STIFF 

500  -  IOOO 

480  •  75  0 

STIFF 

1000  -  2000 

75 0  •  950 

VERY  STIFF 

2000  -  4000 

950  - 1300 

VERT  SOFT 

0  -  2SO 

SOFT 

2S0  -  5 00 

C?!  MW 

STEEL 

MED.  STIFF 

500  -  IOOO 

STIFF 

1000  -  2000 

VERY  STIFF 

2000  -  4000 

PULLOUT 

CAPACITY 

Qu/t  *  CNalTR*  *  CaZTTRI  Tuff  =  CAZTTRL 

Ty/f  U90ER  SUSTAINED  LOAD  MAY  BE 
LIMITED  BY  OTHER  FACTORS. 


ULTIMATE  LOAD  FOR  BUCKLING  of  steel  piles  in  soft  clay. 

Qul*  «  A  (C  E  l)y‘ 

A  =  8  FOR  VERY  SOFT  CLAY .  10  FOR  SOFT  CLAY  I DIMENSIONLESS I 

C  «  SOIL  SHEAR  STRENGTH 

I  -  MOMENT  OF  EMERTIA  OF  CROSS  SECTION 
£  »  MODULUS  OF  ELASTICITY  OF  STEEL. 


}  POUND  AMD 
INCH  UNITS 


FORMULA  APPLIES  TO  SLEHOER  STEEL  PILES.  GENERALLY  HEAVY  STEEL 
PILES  OR  TIMBER  OR  CONCRETE  PILES  ARE  NOT  SUBJECT  TO  BUCKLING 
IF  EMBEDDED  IN  SOIL  FOR  THEIR  ENTIRE  LENGTH. 


Figure  31.  Ultimate  load  capacity  of  piles  in  cohesive  soils  (Ref  11). 
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ANCLE  Of  INTERNAL  FRICTION  BEFORE  PILE  DRIVING. 


DEFINITIONS: 


DEGREES 


PILE  IS  DRIVER  r I  THRU!  JETTl  RG  OR  REMOVAL  OF 
MATERIAL  *1  TRIM  PILE. 

Qyft  a  ULTIMATE  LOAD  CAPACITY.  A  PPLIFD  LOAD  OMLY. 
q>  .  ULTIMATE  BEAR  IMG  CAPACITY  AT  TIP. 

K„*  RATIO  OF  HORIZONTAL  TO  VERTICAL  earth  PRESSURE 
OM  SIDE  OF  PILE  ABOVE  PLASTIC  ZONE.  AVERAGES 
0.5  FOR  f#  •  30°.  1.0  FOR  45°. 

A 'g*RATIO  OF  HORIZOMTAL  TO  VERTICAL  EARTH  PRESSURE 
OM  SIDE  OF  PILE  MIThIN  PLASTIC  ZONE. 

0, ‘ANGLE  OF  I MTERMA  L  FRICTIOM  IM  COMPACTED  ZOME 
AROUMO  PILE  TIP  (4°  TO  S° LARGER  Tham  4). 


i  LIMIT  OF - '  I 

I  SOIL  COMPACTION/ 
DUE  TO  DRIVIH6 


•AVERAGE  PILE  PERIMETER 


Quit  *  %*s*  (  ^  fr  £  ton  cT) 5  -  ■ 

$  ‘  frff*r  *  XefrL/Uj  ~  frL  ' 

9  ■  isk  C> '^RNf  ~  *e  'fr  «-  **;  - , 

EFFECT  OF  t ATEP 
MATER  AT  GROUMD  SURFACE: 

SUBSTITUTE  fsua  ' OR  fT  IN 
FOR  MATER  AT  INTERMEDIATE  LEVEL . 

BE  TREE  N  LIMITING  CONDITIONS. 


IMA  TER  BELOM  d0 1 

FOR  L  i  20/7 

fT  L  FOR  L  <  20 R 


INTERPOLA  TE 


PULLOUT  CAPACITY 


sC&frSton  O') 


COEFF  1C  IENT  Kg 


AV TRANE  PERIMETER  SUPPORTS) 


1000  3000  30CC  WOO 

AVERAGE  OVERBURDEN  PRESSURE  (P.  ) 

IN prF,  FOP  LENGTH  C 


1000  2000  3000  VOOO 

OVERBURDEN  PRESSURE  AT  PUE  VP  (Pt)  W prf 


WHERE  :  0  =  ANGLE  OF  INTERNAL  FRICTION  OF  SO/L 

/C,  •  RATIO  OF  HORIZONTAL  TO  VERTICAL  EARTH  PRESSURE  ON  SIDE  OF  PILE 
TAKEN  TO  BE  EQUAL  TO  I O  IN  GRAPH  (a) 

l'- length  of  firm  embedment 

/Vj  -  VALUES  USED  IN  CHART  @  ARE  TAKEN  FROM  FIGURE  3N 


'WiMr{ 


\  \  \  \n 


w>u«^v|(-;'  t  /Pt 


example  SECT  to  a 


DESIGN  LOAD 


DEFINITIONS 

Qu/f  -  ULTIMATE  PILE  LOAD  CAPACITY. 

-  PERIMETER  SUPPORT  F  END  SUPPORT. 

-  S  X  SURFACE  AREA  *  g  X  END  AREA 

' S‘  VALUES  IN  GRAPH  (?)  APE  FOR  STRAIGHT  CONCRETE  PILES. 

FDR  TAPERED  PILES.  DUE  TO  INCREASED  KH  ,  USE  2*  VALUES  IN  (?)  . 
FOR  STEEL  OP  TIMBER,  OUE  TO  LOWER  COEFFICIENT  OF  FRICTION, 

USE  0.75  K  VALUES  IN  0  . 

r,  |  r,  ,  r,  -  UNIT  WEIGHTS  OF  SOIL,  BY  STRATA. 

(SUBMERGED  vets  below  water  table (wr).) 


WITH  LENGTH  (HI  -  HO'i  L'  -25'j  /}  -  is  pcF  0  -  30- 
AVERAGE  DIAMETER  -  to'  P)-25 pcF  FOR  13 ' 

TIP  AREA  -  0  3HPJ  F*  P?  = BO  pcF  FOR  <?' 

TIMBER  PILE  (CONSIDERED  TAPERED) 

Gulf  «  PERIMETER  SUPPORT  *  END  SUPPORT 
PERIMETER  SUPPORT 

AVERAGE  OVERBUROCN  PRESSURE  (Pa)-(2P;  *  l3/))*(2f,  +  !3r,*2Srs) 

Pc-  129TpvP  ENTER  0  AT  P,  ON  0-30’ 

S  m  ISO  X  2  *  0  7S  »  II2S  puP 
(TAPER)  (TIMBER I 

S  X  SURFACE  AREA  -  II2S  psP  X  TTDL  '  -  73,000  LBS 


FOR  CAPACITY  AS  A 
TENSION  PtU,  USE 
STRAIGHT  PILES  ANO 
USE  PERIMETER  SUPPORT 
ONLY  IN  FARM  EM  BE  DARE NT. 


S  X  SURFACE  AREA  -  U2S  psP  r  jrDL '  *=  73,000  LBS 
END  SUPPORT  - 

OVERBURDEN  PRESSURE  AT  PILE  TIP  (Pt) 

Pt  -  2(7+1317,  *23/)  -  /GO  *  32S  *  IG2S  -  2H0 ppF 
ENTER  0  AT  Pt  •2H0  ANO  0-30’ 

THEN  a  - 39,000  psC 

< %xAs-  38.0 OOpxf  *  OSH  oqPt  -  12,700  LBS 
\QuH  w  73,000  IBS  *  12,  RPOO  IBS  -  86  A  1 


Figure  33.  Ultimate  load  capacity  of  driven  piles  in  cohesionless 
soils  (alternative  method  to  Figure  32)  (Ref  11). 
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Figure  35.  Load  transfer  to  shaft  segment 
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(AL)  versus  downward  movement 
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of  segment. 
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A  contact  element  is  defined  by  two  nodes:  a  spring  constant  or 
"penalty  parameter,"  k,  and  a  fixed  direction  vector,  n.  The  present 
location  of  node  A  (a  =  1,  2)  is  given  by  +  dft,  where  is  the 
initial  position  vector  and  d^  is  the  displacement  vector.  The  contact 
plane  passes  through  the  point  x^  +  d^  and  is  perpendicular  to  n 
(Figure  36a).  The  contact/release  condition  is  defined  as  follows: 

o  >  0 

o  i  0 
where:  a 
£ 

The  quantity  a  is  a  measure  of  the  distance  between  Xg  +  dg  and  the 
contact  plane.  When  contact  is  noted,  a  contact  element  stiffness  and 
out-of-balance  force  are  added  to  the  global  equations. 

If  k  >  0  is  sufficiently  large,  the  point  Xg  +  dg  will  be  forced  to 
lie  (approximately)  on  the  contact  plane.  In  subsequent  steps,  only  the 
stiffness  is  assembled,  and  the  decision  to  remain  in  contact  or  to 
release  is  made  on  the  basis  of  the  sign  of  a,  as  above. 

For  interpreting  output,  the  contact  element  "displacement"  is 
defined  as  o,  and  the  "force"  is  given  by: 

ko  if  o  <  0 

0  if  a  ^  0 

In  addition  to  the  contact  element,  Prevost  has  defined  a  slide 
element.  The  slide-line  element  may  be  used  to  impose  inequality  con¬ 
straints  between  nodes.  Either  perfect  friction  (i.e.,  stick)  or  fric- 
ti onl f-s  (i.e.,  slip)  conditions  may  be  achieved. 

A  slide-line  element  is  defined  by  three  nodes  and  a  spring  con¬ 
stant  or  penalty  parameter,  k.  The  connection  from  node  A  to  node  B 
defines  the  slide-line  direction,  and  node  C  is  the  contact  node 
(Figure  36b). 


(a)  Contact  element. 


n 


(b)  Slide  element. 


Figure  36.  Program  OYNAFLOW. 


The  projected  distance  of  node  C  to  node  A  onto  the  slide- line 
direction  is  denoted  by  a  and  is  given  by: 

a  =  i B  •  fc  /  AB2  0<«H 

where  •  is  the  the  dot  product  of  two  vectors.  The  direction  of  the 
unit  vector  n  to  the  slide-line  direction  is  given  by: 

(1)  in  two  dimensions  by  rotating  A§  as  an  angle  +  n/2 

(2)  in  three-dimensions  by 

n  =  (/te  x  tz)  x  /  (%B  x  tz)  x  fo 

where  x  denotes  the  cross  product  of  two  vectors.  The  local  contact 
stiffness  matrix  k  is  given  by: 

f(l-a)2  a(l  -  a)  -(1  -  a)  ] 

2 

k  =  k  a(l  -  a)  or  -a 

.-(1  -  a)  -a  1  . 

where  the  rows  and  columns  are  arranged  such  that  the  first,  second,  and 
third  rows  (columns)  correspond  to  nodes  A,  B,  and  C,  respectively.  The 
contact/release  condition  is  defined  as  follows:  (1)  in  two-dimensions, 
if  0  %  a  %  1  and  M  0,  then  contact  (otherwise,  release);  and 

(2)  in  three-dimensions,  if  0  £  or  $  1,  then  contact  (otherwise, 
release). 

If  k  is  sufficiently  large,  the  point  C  will  be  forced  to  lie 
(approximately)  on  the  slide-line  AB.  In  subsequent  steps,  only  the 
contact  stiffness  is  assembled,  and  the  decision  to  remain  in  contact  or 
not  is  made  as  described  above. 


Poulos  (Ref  17)  concludes  that  ultimate  load  capacity  and  cyclic 
stiffness  decrease  with  increasing  numbers  of  cycles  and  increasing 
cycle  load  level.  This  becomes  more  significant  when  the  cyclic  load 
approaches  one-half  the  static  ultimate  load.  The  cyclic  degradation 
appears  to  begin  at  the  top  of  the  pile  and  progresses  downward, 
resulting  in  a  gradual  transfer  of  load  to  the  lower  position  of  the 
pile.  The  crucial  factor  in  determining  the  amount  of  cyclic  degrada¬ 
tion  is  the  shear  strain  for  skin  friction. 

Poulos  (Ref  18)  conducted  a  limited  investigation  of  an  effective 
stress  approach  to  determine  pore  pressure  increases  with  cyclic  loading 
and  the  resulting  modulus  degradation  factors.  Possibly,  a  lack  of  data 
prevented  further  development.  Figure  37  shows  a  compilation  of  ob¬ 
servations  (Ref  17)  showing  degradation  as  a  function  of  strain  ratio, 
where: 

yss  =  (0.10  to  0.25)  ys 

and 

Y  =  static  shear  strain  to  failure 
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Figure  37.  Distribution  of  force 
within  pile,  75,  150 
and  300  kg  loads. 
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Poulos  (Ref  16)  reports  that  "one-way  cyclic  loading  (i.e. ,  cyclic 
loading  between  zero  minimum  load  and  a  specified  maximum)  produces 
significantly  less  degradation  of  pile  capacity  than  would  occur  with 
"two-way"  loading  (i.e.,  loading  alternating  between  tension  and  com¬ 
pression  with  zero  as  a  mean  value).  He  points  out  that  degradation 
will  occur  at  different  rates  along  the  pile,  depending  on  local  stress 
level.  Even  in  initially  homogeneous  soil,  a  nonuniform  distribution  of 
soil  modulus  and  skin  friction  will  result  from  cyclic  loading  because 
nonuniformity  of  stress  distribution  occurs  along  the  pile.  Degradation 
occurs  in  the  ultimate  skin  friction  along  the  length  of  the  pile  and 
also  in  the  ultimate  base  resistance.  The  major  problem  in  a  cyclic 
response  is  determining  how  the  degradation  factors  vary  with  strain  and 
number  of  cycles. 

The  cyclic  shear  strain  in  the  soil  adjacent  to  the  pile,  y  ,  can 
be  estimated  as: 


where:  Pc  =  cyclic  displacement  of  pile  at  a  point  on  the  pile  shaft 
d  =  pile  diameter 
l  =  An  [5Y  (1  -  ms)  L/d] 

f  =  1.0  homogeneous  infinitely  deep  soil  and 
0.5  modulus  increases  with  depth 

Ps  =  Poisson's  ratio 

L  =  embedded  pile  length 


The  cyclic  base  strain  can  be  estimated  from  the  above  assuming 
Ps  =  0.5,  L/d  =  100. 


0  •  4  Pt 


Yc  = 


where  P^  is  the  cyclic  displacement  of  pile  base  and  d^  is  the  diameter 
of  pile  base. 


It  is  important  to  note  that  the  above  is  based  on  the  assumptions 
of  elastic  behavior.  Typical  predictions  for  pile  settlement  take  the 
form  of  Figure  38.  Also  shown  in  Figure  38  is  the  ultimate  cyclic  load 
as  a  function  of  number  of  cycles.  The  settlement  problem  is  the  major 
concern  for  piles  in  cohesionless  materials. 

Friction  Between  Calcareous  Sand  and  Building  Materials 

Calcareous  sediments  have  proven  troublesome  to  offshore  facil¬ 
ities.  Piles  driven  into  calcareous  sands  have  been  noted  to  penetrate 
and  to  be  extracted  with  much  less  effort  than  predicted  by  conventional 
techniques.  A  research  program  was  previously  conducted  at  NCEL 
(Ref  19)  in  which  calcareous  sediments  were  collected  from  three  envi¬ 
ronments:  a  deep-ocean  site,  a  shallow-ocean  site,  and  coralline  sand 
from  an  atoll  beach.  The  coefficients  of  friction  of  these  sands  and  of 
a  quartz  sand  used  as  a  standard  were  measured  against  surfaces  of  rough 
and  smooth  steel  and  mortar.  Volume  changes  were  measured  as  a  function 
of  sliding  displacement. 

Experience  has  caused  engineers  to  reduce  pile  capacities  in  cal¬ 
careous  materials.  This  usually  results  in  load  capacity  reductions  to 
one-fourth  that  of  piles  in  normal  materials.  In  the  past,  it  was  not 
clearly  understood  why  calcareous  materials  exhibit  inferior  pile  sup¬ 
port.  Measured  angles  of  internal  friction  --  a  measure  of 
strength  —  are  high:  34  degrees  or  greater.  Part  of  the  problem  was 
identified  as  the  low  increase  in  soil  effective  stress  during  pile 
driving,  which  is  thought  to  result  from  a  crushing  or  collapse  of  a 
cemented  soil  structure  or  from  the  breakup  of  individual  carbonate 
grains.  This  lack  of  increase  in  effective  stress  results  in  a  rela¬ 
tively  lower  shear  strength  in  the  soil  mass  surrounding  the  pile. 
Further,  it  was  thought  that  the  coefficients  of  friction  might  be 
substantially  less  between  calcareous  materials  and  piles.  Valent 
(Ref  18)  conducted  a  test  program  using  a  modified  direct  shear  test 
machine  in  which  the  lower  half  of  the  apparatus  contained  the  building 
material  and  the  upper  half  the  soil  sample.  Table  2  summarizes  his 
results. 
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(c)  Number  of  cycles. 

Figure  38.  Pile  degradation  with  cyclic  loading  (Ref  18) 


Table  2.  Summary  of  Friction  Test  Results  (Ref  18) 


Test  No. 

Soil  Material 

a 

^peak 

a 

^residual 

Base  Material*3: 

Sandc 

1 

Quartz  sand 

0.  67d 

0.54 

2 

Coralline  sand 

0.66 

0.56 

3 

Coralline  sand 

0.68. 

0.57 

4 

Oolitic  sand 

0.77d 

0.61 

5 

Oolitic  sand 

0.81 

0.62 

6 

Foram  sand-silt 

0.  64d 

0.58 

Base  Material*3:  Smooth  Steel 

7 

Quartz  sand 

0.  27d 

0.19 

8 

Coralline  sand 

0.20 

°-  !7 

9 

Coralline  sand 

0.20 

0. 18e 

10 

Coralline  sand 

0.21. 

0.17 

11 

Oolitic  sand 

0. 15d 

0.13 

12 

Oolitic  sand 

0.32 

0.31 

13 

Foram  sand-silt 

0.40 

0.37 

Base  Material*3:  Rough  Steel 

14 

Quartz  sand 

0.60 

0.54 

15 

Coralline  sand 

0.63 

0.55 

16 

Oolitic  sand 

0.54 

0.51 

17 

Oolitic  sand 

0.58, 

0.50 

18 

Foram  sand-silt 

0.66 

Base  Material*5:  Smooth  Concrete  . 

19 

Quartz  sand 

0.60 

0.54 

20 

Coral  1 i ne  sand 

0.63 

0.56 

21 

Oolitic  sand 

0.59 

0.52 

22 

Oolitic  sand 

0.58, 

0.54 

23 

Foram  sand-silt 

0.67 

Base  Material*3:  Rough  Concrete 

24 

Quartz  sand 

0.69 

0.57 

25 

Coralline  sand 

0.66 

0.59 

26 

Oolitic  sand 

0.74 

0.57 

aFor  direct  shear  tests,  y  =  tan  $  where  <)>  =  angle  of  internal 
friction;  for  friction  tests,  y  =  tan  6  where  6  =  angle  of 
sliding  friction. 

^Soil  in  bottom  shear  ring  for  direct  shear  tests,  or  building 
material  in  friction  tests. 

cBase  material  same  as  soil  material  for  direct  shear  tests. 

aThese  tests  run  with  mechanical  measurement  system; 
i.e.,  proving  ring  and  manual  recording  of  data. 

etow  value  for  y  reached  shortly  after  y  .  ,  thereafter  y 
increased  with  displacement  to  end  of  t£f?K 

fNo  peak  y  reached,  y  increasing  through  end  of  test. 


In  general,  the  results  show  that  the  low  friction  forces  in  cal¬ 
careous  sediments  are  not  the  result  of  low  achievable  coefficients  of 
friction  between  calcareous  sediments.  Since  the  coefficients  of  fric¬ 
tion  of  calcareous  sands  are  comparable  to  other  sands,  low  friction, 
then,  must  be  attributable  to  low  normal  force.  Deep-ocean  sand  (foram- 
iniferal  sand-silt)  exhibits  one  possible  cause  for  low  developed  normal 
force.  The  volume  change  during  testing  indicated  a  considerable  volume 
decrease  during  development  of  resisting  friction  force,  probably  due  to 
crushing  of  the  skeletal  structures  and  shell  fragments.  Penetration  of 
a  pile  in  such  a  material  would  crush  the  hollow  shell  material  with 
only  a  minimal  increase  in  effective  stress  of  the  surrounding  material. 

Quoting  from  Valent  (Ref  19): 


1.  The  calcareous  sediments  tested,  and  presumably  calcareous 
sediments  in  general,  develop  coefficients  of  friction  against 
steel  and  concrete  building  materials  that  are  comparable  to 
those  developed  by  quartz-type  sands.  Thus,  the  possibility 
of  low  coefficients  of  friction  being  responsible  for  the 
observed  low  friction  forces  on  driven  piling  and  other  pene- 
trators  in  calcareous  materials  is  ruled  out. 

2.  The  observed  large  volume  decreases  during  shear  of  the 
foraminiferal  sand-silt  are  probably  responsible  for  the  low 
developed  friction  forces  in  these  hollow-shelled  materials. 
Such  large  volume  decreases  at  nonincreasing  normal  load  imply 
densif ication  in  the  field  without  accompanying  increases  in 
normal  stress  on  the  penetrator  surface. 

3.  Low  developed  friction  forces  in  other  calcareous  mate¬ 
rials  may  arise  from  a  similar  mechanism  involving  a  hypoth¬ 
esized  loose,  but  cemented,  structure  for  the  soil  material. 
The  application  of  shear  stresses  during  penetration  would 
cause  collapse  of  this  structure  to  a  denser,  but  still  loose, 
arrangement. 


Piles  in  Calcareous  Sands 


Calcareous  sands,  as  discussed,  are  noted  for  loose  arrangement  of 
particles  lightly  cemented  to  form  a  structure  to  support  other  layers 
without  compacting.  However,  upon  shearing,  the  structure  is  destroyed, 


breaking  the  cement  bonds;  the  loose-grained  structure  then  compacts  and 
densities.  The  sand  (Lade  sand)  is  a  very  loose  sand  (30%  relative 
density)  and,  as  such,  is  a  "manufactured"  sand  that  would  not  occur  in 
nature.  The  properties  of  this  material  might  qualitatively  be  expected 
to  represent  the  constitutive  behavior  of  a  calcareous  sand. 

To  analyze  a  pile's  load  capacity  using  the  effective  stress  soil 
model,  the  mesh  in  Figure  39  was  used.  This  mesh  simulates  a  pile 
already  in  place  (i.e.,  not  the  driving  of  the  pile).  Horizontal  and 
vertical  springs  were  used  to  join  the  pile  to  the  soil  field  to  simu¬ 
late  the  interface  and  allow  for  pile  movement.  Unfortunately,  DYNAFLOW 
permits  only  elastic  springs  at  this  time. 

The  soil  properties  used  wero  those  of  Lad  sand.  The  material 
properties  simulate  the  condition  of  the  soil  after  placement  of  the 
pile.  No  densification  occurs,  since  the  pile  is  in  place  at  the  start 
of  the  analysis.  A  prediction  of  pile  capacity  was  made  by  using  con¬ 
ventional  techniques  found  in  Figures  32  and  33  adjusted  to  the  problem 
conditions  but  without  densification.  Results  indicated  a  pile  resis¬ 
tance  in  normal  average  cohesionless  sand  of  at  least  1,500  kg  with 
about  150  kg  of  that  in  side  friction  and  the  remainder  in  end  bearing. 

Load  Settlement  of  the  Pile.  Figure  40  shows  the  load  settlement 
of  the  pile.  Several  values  we»e  tried  for  the  spring  constants. 
Results  show  the  pile  experiences  large  settlements  between  300  and 
400  kg.  Figure  37  shows  the  distribution  of  force  within  the  pile  for 
three  load  levels.  Note  that  most  of  the  increase  in  load  is  due  to 
skin  friction.  Figure  41  shows  the  distribution  of  force  along  the  pile 
at  a  point  in  the  loading  having  a  force  of  300  kg  applied  at  the  top  of 
the  pile.  The  results  in  Figure  41  show  about  one-half  to  two-thirds  of 
the  load  is  taken  in  skin  friction.  The  skin  friction  portion  of  the 
load  agrees  with  the  predicted  values;  however,  the  end  bearing  is 
substantially  less.  It  is  important  to  note  that  the  sand  used  (Lade 
sand)  had  a  relative  density  of  only  30%  which  is  so  loose  it  does  not 
occur  in  nature  but  is  "manufactured"  in  the  laboratory.  The  level  of 
loading  is  about  one-fourth  that  of  a  pile  in  normal  sand,  which  is  the 
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level  expected  from  experience  in  calcareous  sands.  Figure  42  shows 
contours  of  principal  stress  around  the  pile  at  a  typical  load  step. 
These  show  typical  patterns,  as  expected;  the  level  of  stress  below  the 
pile  tip  would  be  on  the  last  surface  at  yield. 

Results  show  that  the  order  of  magnitude  of  the  pile  capacity  in 
calcareous  materials  is  predicted  correctly.  Further,  the  friction 
developed  on  the  sides  of  the  pile  is  at  the  expected  level,  not  sub¬ 
stantially  reduced  from  normal  sands.  The  amount  of  end  bearing  is 
slightly  dependent  on  the  vertical  spring  constant.  Use  of  nonlinear 
springs  would  have  allowed  slip  to  occur,  transferring  more  load  into 
end  bearing  when  exceeding  some  local  slip  level. 

Pile  Behavior  Under  Cyclic  Loading.  Next,  pile  behavior  under 
cyclic  loading  was  examined.  The  applied  loading  consisted  of  cyclic 
variation  of  the  vertical  load  on  the  pile.  In  one  case,  the  loading 
was  cycled  in  "one-way"  loading  (zero  minimum  and  specified  maximum)  and 
in  "two-way"  loading  (alternating  compression  and  tension).  The  first 
example  was  at  a  relatively  low  level  of  loading  for  the  model  pile 
(about  one-seventh  of  yield).  Figure  43  shows  the  stress  beneath  the 
pile  tip  for  one-way  loading  (0  to  50  kg);  Figure  44  shows  the  shear 
stress  contours  around  the  pile.  Figure  45  shows  the  stress  beneath  the 
pile  tip  for  two-way  loading  (-50  to  50  kg);  the  shear  stress  plot  is 
similar  to  that  of  Figure  46.  Note  that  the  tip  forces  reduce  with 
loading,  transferring  more  to  friction  to  maintain  load  levels. 
Figure  46  shows  the  typical  contact  element  force,  a  measure  of  friction 
between  the  soil  and  pile.  Note  the  gradual  buildup  in  the  two-way 
loading  case. 
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Figure  43.  Horizontal  and  vertical  stress  of  pile  in  soil 
field  -  cyclic  load,  0  to  50  kg  (soil  element 
beneath  pile  tip). 
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Figure  46.  Friction  force  in  pile  in  soil  field  with  cycle  loads  of 
0  to  50  and  -50  to  50  kg. 


The  loading  was  increased  to  about  40%  of  yield.  Figure  47  shows 
the  distribution  of  force  within  the  pile  for  both  one-way  and  two-way 
loading  for  the  first  and  fifth  cycles.  Note  the  increase  in  friction 
loading  in  the  pile  with  two-way  loading,  showing  the  degradation  with 
cyclic  loading.  Note  also  that  for  both  cases  the  friction  increases 
more  in  the  lower  half  of  the  pile.  Figure  48  composes  the  typical 
contact  element  force,  again  showing  the  increase  in  friction  with  each 
cycle.  Figures  49  and  50  show  the  stress  beneath  the  pile  tip.  Both 
drop  off  with  each  cycle;  however,  the  two-way  loading  does  so  at  a 
faster  rate.  Figure  51  shows  the  pile  settlements  for  both  cases.  Note 
that  the  change  in  settlement  between  the  first  and  fifth  cycle  is  over 
four  times  greater  for  the  two-way  loading.  Figure  5?  shows  the  de¬ 
formed  mesh  for  the  one-way  loading;  the  mesh  would  be  similar  for  the 
two-way  loading. 

The  loading  was  increased  to  80%  of  pile  capacity,  and  the  one-way 
loading  case  was  tried.  This  load  exceeds  the  friction  capacity  by 
itself,  so  two-way  loading  could  not  be  used.  Figure  53  shows  the 
stress  beneath  the  pile  tip,  and  Figure  54  shows  the  deformed  mesh. 
Failure  occurs  between  steps  70  and  80.  Figure  55  shows  the  shear 
stress  contours  around  the  pile  at  steps  10  and  70.  More  stress  is 
transferred  from  tip  end  bearing  to  side  friction. 
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Figure  50.  Pile  in  soil  field  -  cyclic  load,  -150  to  150  kg  (soil 
element  beneath  pile  tip). 
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Figure  51.  Pile  in  soil  field  -  cyclic  load  (pile  settlement). 
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Figure  54.  Deformed  mesh  -  cyclic  load,  0  to  300  kg. 
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Figure  54. 
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The  cyclic  results  clearly  show  the  degradation  of  the  pile  under 
cyclic  loading.  Results  show  the  two-way  loading  causes  more  degrada¬ 
tion,  as  expected.  The  model  appears  to  be  performing  well  in  predic¬ 
ting  the  qualitative  pile  behavior,  particularly  the  settlements 
(Figures  37  and  51). 

The  pile  problems  were  repeated  for  the  undrained  case  in  which  the 
water  table  was  at  the  surface.  Figure  56  shows  the  load  settlement 
curve.  The  original  load  increment  used  in  the  drained  problem  above 
was  too  coarse  for  the  undrained  problem  and  produced  an  instability  at 
yielding.  The  load  step  was  reduced,  and  the  solution  proceeded  in  a 
satisfactory  manner.  The  pile  capacity  was  about  250  kg,  lower  than  the 
350  kg  in  the  drained  case.  The  loading  was  cycled  in  one-way  and 
two-way  conditions.  Figure  57  presents  the  distribution  of  force  within 
the  pile,  showing  the  increase  in  friction  with  cycling,  as  with  the 
drained  case.  Friction  is  slightly  greater  in  this  case.  Figures  58 
and  59  show  the  soil  stress  beneath  the  pile.  The  end  bearing  stress 
drops  off  with  cycling,  showing  the  shift  to  friction.  The  two-way 
loading  has  a  faster  degradation,  as  clearly  shown  in  Figure  60,  which 
shows  the  settlement.  The  two-way  loading  produces  a  more  permanent 
settlement  differential  after  the  fifth  cycle.  Figure  61  shows  the  pore 
pressure  in  the  soil  beneath  the  pile  tip.  The  pore  pressure  rises  to 
about  three-quarters  of  the  continuing  stress  and  then  drops  off  as  the 
load  is  shifted  to  friction  on  the  sides.  Figure  62  shows  the  friction 
force.  Note  the  degradation  effects  with  cycling  as  load  drops  off  with 
each  cycle. 
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Figure  58.  Pile  in  soil  field  *  cyclic  load,  0  to  150  kg  (soil  element 
beneath  pile  tip)  (with  water  table). 
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Figure  69.  Pile  in  soil  field  -  cyclic  load,  -160  to  150  kg  (soil 
element  beneath  pile  tip)  (with  water  table). 
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Figure  60.  Pile  in  soil  field  -  cyclic  load  (pile  settlement)  (with 
water  table). 
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Figure  61.  Pile  in  soil  field  -  cyclic  load  (soil  element  beneath  pile 
tip)  (with  water  table). 
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Figure  62. 


-  Cyclic  load  0  to  ISO  kg 

- Cyclic  load,  -150  to  150  kg 


Pile  in  soil  field  -  friction  force  in  pile  (with  water 
table). 


CONCLUSION 

In  this  work  the  accuracy  and  ability  of  the  Prevost  soil  model  has 
been  compared  to  laboratory  test  data.  Results  show  favorable  agree¬ 
ment.  The  model  has  been  Implemented  Into  a  finite  element  code  of 
considerable  significance.  The  finite  element  formulation  allows  for 
computation  of  effective  stresses,  pore  pressures,  and  fluid  flow  under 
static  and  transient  loadings.  The  material  properties  required  to 
utilize  the  Prevost  soil  model  are  determined  by  conventional  static 
triaxial  compression  and  extension  tests.  Once  these  test  data  have 
been  obtained,  the  formulations  of  required  input  parameters  to  quantify 
the  properties  of  the  soil  are  determined  by  an  automated  procedure. 

Several  demonstration  problems  Illustrating  the  cyclic  degradation 
capability  of  the  soil  model  have  been  performed  with  successful  re¬ 
sults.  At  this  point  the  model  has  been  shown  to  agree  reasonably  well 
with  laboratory  tests  which  impose  known  simple  loading  states.  Actual 
data  of  field  behavior  is  very  limited.  The  demonstration  cases  show 
the  model  performs  consistently  with  field  experience;  however,  more 
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quantitative  data  coaiparison  is  required.  The  soil  properties  are 
dependent  upon  the  stress- state  existing,  which  is  a  function  of  the 
type  of  loading  condition  or  structure  constructed  at  a  site.  Thus, 
actual  case  studies  must  have  the  complexity  of  evaluating  material 
properties  correctly  before  evaluation  of  response  under  load  can  be 
evaluated.  In  addition,  verification  of  the  model  under  more  complex 
stress  states  in  which  drainage  occurs  is  required. 

It  is  recommended  that  this  procedure  be  tried  in  a  field  study 
where  actual  field  data  can  be  obtained  and  utilized  in  a  demonstration 
boundary  value  case  study.  Validation  of  the  model  and  code  in  a  quan¬ 
tified  environment  should  be  performed.  A  suitable  bench  mark  must  be 
obtained  upon  which  analytical  results  can  be  compared. 
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(LAWRENCE  LIVERMORE  LAB.  TOKARZ);  La  Jolla  CA  (Acq  Dept.  Lib.  C-U75AI:  M  Duncan. 
Berkeley  CA 

UNIVERSITY  OF  DELAWARE  Newark.  DE  (Dept  of  Civil  Engineering.  Chesson) 

UNIVERSITY  OF  HAWAII  HONOLULU.  HI  (SCIENC  E  AND  LECH  DIV.) 

UNIVERSITY  OF  ILLINOIS  (Hall)  Urbana.  IL;  Metz  Ref  Rm.  Urbana  IL.  URBANA.  11  (DAVISSON); 

URBANA.  IL  (LIBRARY);  Urbana  IL  (CE  Dept.  W.  Gamble) 

UNIVERSITY  OF  MASSACHUSETTS  (Heroncmus).  ME  Dept.  Amherst.  MA 
UNIVERSITY  OF  MICHIGAN  Ann  Arbor  Ml  (Riehart) 

UNIVERSITY  OF  NOTRE  DAME  Katona.  Notre  Dame.  IN 
UNIVERSITY  OF  TEXAS  Inst.  Marine  Sei  (Library).  Port  Arkansas  TX 

UNIVERSITY  OF  TEXAS  AT  AUSTIN  AUSTIN.  IX  (THOMPSON):  Austin  TX  (R.  Olson):  Austin.  TX 
(Breen) 

UNIVERSITY  OF  WASHINGTON  Dept  of  Civil  Engr  (Dr.  Mattoek).  Seattle  WA.  SEATTLE.  W'A  (OCEAN 
ENG  RSCH  LAB.  GRAY);  Seattle.  WA  Transportation.  Construction  &  Geom.  Div 
AMETEK  Offshore  Res.  &  Engr  Div 
ATLANTIC  RICHFIELD  CO.  DALLAS.  TX  (SMITH) 

BECHTEL  CORP  SAN  FRANCISCO.  CA  (PHELPS) 

BETHLEHEM  STEEL  CO.  Dismuke.  Bethelehem.  PA 

BRITISH  EMBASSY  M  A  Wilkins  (Sei  &  Tech  Dept)  Washington.  DC 

BROWN  &  ROOT  Houston  TX  (D  Ward) 

CHAS.  TL  MAIN.  INC.  (R.C.  Goveltc).  Portland.  OR 

CHEVRON  OIL  FIELD  RESEARCH  CO.  LA  HABRA.  CA  (BROOKS) 

CONRAD  ASSOC.  Van  Nuys  CA  (A  Luisoni) 

DR  A  VO  CORP  Pittsburgh  PA  (Wright) 

EVALUATION  ASSOC  INC  KING  OF  PRUSSIA.  PA  (FEDELE) 

EXXON  PRODUCTION  RESEARCH  CO  Houston.  TX  (Chao) 

FURGO  INC.  Library.  Houston.  TX 

GEOTECHNICAL  ENGINEERS  INC.  Winchester.  MA  (Paulding) 

GRUMMAN  AEROSPACE  CORP  Bcthpage  NY  (Tech  Info.  Ctr) 

HALEY  &  ALDRICH.  INC  Cambridge  MA  (Aldrich.  Jr.) 

NUSC  DET  Library.  Newport.  Rl 

LAMONT-DOHERTY  GEOLOGICAL  OBSERVATORY  Palisades  NY  (McCoy) 

MARATHON  OIL  CO  Houston  TX 

MOBIL  R  &  D  CORP  Manager.  Offshore  Engineering,  Dallas.  TX 

MUESER.  RUTLEDGE.  WENTWORTH  AND  JOHNSTON  New  York  (Richards) 

PORTLAND  CEMENT  ASSOC.  Skokie  IL  (Rsch  &  Dev  Lab.  Lib.) 

RAND  CORP.  Santa  Monica  CA  (A.  Laupa) 

RAYMOND  INTERNATIONAL  INC.  E  Colic  Soil  Tech  Dept.  Pennsauken.  NJ 
SANDIA  LABORATORIES  Library  Div..  Livermore  CA 
SCHUPACK  ASSOC  SO  NORWALK.  CT  (SCHUPACK) 

SHANNON  &  WILLSON  INC.  Librarian  Seattle.  WA 
SHELL  DEVELOPMENT  CO.  Houston  TX  (E.  Doyle) 

SHELL  OIL  CO  HOUSTON.  TX  (MARSHALL) 

TRW  SYSTEMS  CLEVELAND.  OH  (ENG  LIB  ):  REDONDO  BEACH.  CA  (DAI) 
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DEPARTMENT  OF  THE  NAVY 

NAVAL  Civil  ENGINEERING  LABORATORY 
port  hueneme.  California  93043-5003 

OFFICIAL  BUSINESS 

PENALTY  FOR  PRIVATE  USE.  SSOO 
I  IMO-NCEL  2700/4  (REV  12-73) 
OS30-LL-L70-0044 


Commanding  Officer 
Code  LI  4 

Naval  Civil  Engineering  Laboratory 
Port  Hueneme,  California  93043-5003 


INSTRUCTIONS 


The  Naval  Civil  Engineering  Laboratory  has  revised  its  primary  distribution  lists.  The  bottom  of 
the  mailing  label  has  several  numbers  listed.  These  numbers  correspond  to  numbers  assigned  to  the  list  of 
Subject  Categories.  Numbers  on  the  label  corresponding  to  those  on  the  list  indicate  the  subject  category  and 
type  of  documents  you  are  presently  receiving.  If  you  are  satisfied,  throw  this  card  away  (or  file  it  for  later 
reference). 

If  you  want  to  change  what  you  are  presently  receiving: 

•  Delete  -  mark  off  number  on  bottom  of  label. 

•  Add  -  circle  number  on  list. 

•  Remove  my  name  from  all  your  lists  -  check  box  on  list. 

•  Change  my  address  -  line  out  incorrect  line  and  write  in  correction  (ATTACH  MAILING  LABEL). 

•  Number  ot  copies  should  be  entered  after  the  title  of  the  subject  categories  vou  select 
fold  on  line  below  and  drop  in  the  mail. 

Note  Numbers  on  label  but  not  listed  on  questionnaire  are  for  NCEL  use  only,  please  ignore  them. 
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NAVAL  CIVIL  ENGINEERING  LABORATORY 
PORT  HUENEME.  CALIFORNIA  93043 
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PENALTY  POP  private  USE.  $300 
1  IND.NCEL.2700/ 4  (REV.  12-73) 

OS  3Q.LL-L7Q.Q044 


Commanding  Officer 
Code  LI  4 

Naval  Civil  Engineering  Laboratory 
Port  Hueneme,  California  93043 


DISTRIBUTION  QUESTIONNAIRE 

The  Naval  Civil  Engineering  Laboratory  is  revising  its  primary  distribution  lists. 


SUBJECT  CATEGORIES 

1  SHORE  FACILITIES 

2  Construction  methods  and  materials  (including  corrosion 

control,  coatings) 

3  Waterfront  structures  (maintenance/deterioration  control) 

4  Utilities  (including  power  conditioning) 

5  Explosives  safety 

6  Construction  equipment  and  machinery 

7  Fire  prevention  and  control 

8  Antenna  technology 

9  Structural  analysis  and  design  (including  numeric* I  and 

computer  techniques) 

10  Protective  construction  (including  hardened  shelters, 

shock  and  vibration  studies) 

1 1  Soil /rock  mechanics 

13  BEQ 

14  Airfitlds  and  pavements 

15  ADVANCED  BASE  AND  AMPHIBIOUS  FACILITIES 

16  Base  facilities  (including  shelters,  power  generation,  water  supplies) 

17  Expedient  roads/air  fields/bridges 

18  Amphibious  operations  (including  breakwaters,  wave  forces) 

19  Over-thc-Beech  operations  (including  containerisation, 

materiel  transfer,  lighterage  and  cranes) 

20  POL  storage,  transfer  and  distribution 
24  POLAR  ENGINEERING 

24  Same  as  Advanced  Base  and  Amphibious  Facilities, 
except  limited  to  cold-region  environments 


TYPES  OF  DOCUMENTS 

85  Techdata  Sheets  86  Technical  Reports  and  Technical 

83  Table  of  Contents  8c  Index  to  TDS 


28  ENERGY/POWER  GENERATION 

29  Thermal  conservation  (thermal  engineering  of  buildings,  HVAC 

systems,  energy  loss  measurement,  povwr  generation) 

30  Controls  and  alectricel  conservation  (electrical  systems. 

energy  monitoring  and  control  systems) 

31  Fuel  flexibility  (liquid  fueit,  coal  utilization,  energy 

from  solid  waste) 

32  Alternate  energy  source  (geothermal  power,  photovoltaic 

power  systems,  solar  systems,  wind  systems,  energy  storage 
systems) 

33  Site  data  and  systems  integration  (energy  resource  data,  energy 

consumption  data,  integrating  energy  systems) 

34  ENVIRONMENTAL  PROTECTION 

35  Solid  waste  management 

36  Hazardous/toxic  materials  management 

37  Wastewater  management  and  sanitary  engineering 

38  Oil  pollution  removal  and  recovery 

39  Air  pollution 

40  Noise  abatement 

44  OCEAN  ENGINEERING 

45  Seafloor  soils  and  foundations 

46  Seafloor  construction  systems  and  operations  (including 

diver  and  manipulator  tools) 

47  Undersea  structures  and  materials 

48  Anchors  and  moorings 

49  Undersea  power  systems,  electromechanical  cables. 

and  connectors 

50  Pressure  vessel  facilities 

51  Physical  environment  (including  site  surveying) 

52  Ocean-based  concrete  structures 

53  Hyperbaric  chambers 

54  Undersea  cable  dynamics 
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